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Abstract 

We compute the effect of soft-gluon resummation, at the next-to-leading-logarithmic level, 
in the hadroproduction cross-section for heavy flavours. Applications to top, bottom and 
charm total cross-sections are discussed. We find in general that the corrections to the fixed 
next-to-leading-order results are larger for larger renormalization scales, and small, or even 
negative, for smaller scales. This leads to a significant reduction of the scale-dependence of 
the results, for most experimental configurations of interest. 
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1 Introduction 



Large logarithmically-enhanced corrections due to soft-gluon radiation are a general feature in 
the study of the production cross-sections of high-mass systems near threshold. Techniques for 
resumming these corrections have been developed over the past several years, starting from the 
case of Drell-Yan (DY) pair production Jl|, |2|. In this paper we deal with the resummation of 
logarithmically-enhanced effects in the vicinity of the threshold region in the total cross-section 
for the hadropro duct ion of heavy quarks. This topic has been examined in several studies in recent 
years 0-0. These analyses were limited to soft-gluon resummation to leading-logarithmic accu- 
racy. The work we present here extends and updates the results of refs. [Q and @ by consistently 
including all the next-to-leading logarithmic contributions. 

The resummation program of soft-gluon contributions is best carried out in the Mellin-transform 
space, or iV-space. In iV-moment space the threshold-production region corresponds to the limit 
iV — > oo and the typical structure of the threshold corrections is as follows 

{oo In ~| 

l + E<E^, m ln m iV , (1) 
n=l m=l ) 

where <xjy°' ) is the corresponding partonic cross-section at leading order (LO). In the case of 
hadronic collisions, the best studied process [[!], [2| is Drell-Yan (DY) lepton-pair production. In 
the DY process the logarithmic terms in the curly bracket of eq. (HD can be explicitly summed and 
organized in a radiative factor Ady,n that has an exponential form: 

{oo n+1 ^ 

E«"E G «- lnmiV ( 2 ) 
n=l m=l ) 

= ex V {\nNg^l(a s \nN)+g^ Y (a s \nN) + a s g§ Y (a s \nN) + ...} . (3) 

Note that the exponentiation in eq. (0) is not trivial. The sum over m in eq. ([T]) extends up to 
m = 2n while in eq. (|2|) the maximum value for m is smaller, m < n + 1. In particular, this 
means that all the double logarithmic (DL) terms OL n s c n ^ n hi 2n iV in eq. ([I]) are taken into account 
by simply exponentiating the lowest-order contribution ol s C\^ In 2 N. Then, the exponentiation in 
eq. (H) allows one to define the improved perturbative expansion in eq. ©. The function In N g^ Y 
resums all the leading logarithmic (LL) contributions a™ ln n+1 N, g^ Y contains the next-to-leading 
logarithmic (NLL) terms a™ ln n N, cn s g DY contains the next-to-next-to-leading logarithmic (NNLL) 
terms ln n N, and so forth. Once the functions g^ Y have been computed, one has a systematic 
perturbative treatment of the region of N in which a s lnN ^1, which is much larger than the 
domain a s In 2 iV 1 in which the fixed-order calculation in a s is reliable. 

The physical bases for the soft-gluon exponentiation in the DY process are dynamics and kine- 
matics factorizations [l], [|. Dynamics factorization follows from gauge invariance and unitarity: 
in the soft limit multi-gluon amplitudes fulfil generalized factorization formulae given in terms 
of a single-gluon emission probability. Kinematics factorization follows from the fact that the 
multi-gluon phase space in the appropriate soft limit can exactly be factorized by working in iV 
space. 

The general extension of these results to the resummation of NLL contributions in heavy- 
quark production was first considered by Kidonakis and Sterman || and more recently examined in 



1 



ref. |9| . There are kinematics and dynamics complications in going from the DY to the heavy-quark 
processes. Although in the heavy-quark case the kinematics is more involved it can nonetheless be 
recast in factorized form in iV-moment space |8|, [|. Dynamics complications are more essential. 
While in the DY case resummation is only associated with the gg-annihilation subprocess, in the 
heavy-quark case one has to consider both the qq and gg subprocesses. Moreover, unlike the DY 
vector boson, the heavy-quark pair carries colour charge. Thus the heavy-quark cross-section is 
sensitive to radiative corrections due to soft-gluon emission from the heavy quark and antiquark 
produced in the final state. In particular, this emission induces colour exchange in the hard 
subprocess, whose colour state can thus fluctuate. The presence of colour fluctuations implies 
that, after soft-gluon resummation, the partonic cross-sections have an exponentiated form that 
is more involved than the simple exponential in eq. (J3|). 

The general structure of the resummed partonic cross-section for heavy-quark hadroproduction 
is as follows || 

=E M Ln Mw, N , (4) 

i,j 

where all the various factors depend on the heavy-quark mass m and the renormalized coupling 
a s (fi 2 ) evaluated at the scale /i, although this dependence is not explicitly denoted. The sum in 
eq. (f|) is over all possible colour states I, J of the LO hard-scattering subprocess i+j — > QQ. The 
colour amplitude My^jv is computable as a power series expansion in a s , but it does not contain 
any logarithmically-enhanced corrections produced by soft-gluon radiation. These corrections are 
embodied by the factor A^jv that is a matrix in the space of the colour states. Performing the all- 
order resummation, the factor A^jv can be expressed in exponential form but, unlike A dy ,n, we 
are now dealing with an exponential matrix and exponentiation has to be understood in a formal 
sense. The analogues of the functions g^ Y are colour matrices. They have to be evaluated and 
then diagonalized, before one can rewrite the right-hand side of eq. ([|) as an actual exponential. 

The resummation formula (|4]) applies to the general case of heavy-quark production at fixed 
invariant mass of the produced quark pair and at fixed scattering angle In this case, the radiative 
factor Ajj ^v depends on the difference in rapidity between the heavy quark and antiquark and the 
NLL resummation has been worked out explicitly in ref. |§. 

The goal of this paper is the resummation of logarithmically-enhanced corrections to the heavy- 
quark total hadroproduction cross-section, up to NLL order. The non-trivial angular dependence 
of the eigenvalues and eigenvectors of the radiative factor A^jv calculated in ref. ||, however, 



makes it difficult to apply in practice [11] this approach to the evaluation of the total production 
cross-section. A new general method for the resummation of soft-gluon contributions to hard- 
scattering processes has recently been developed in refs. |], This technique allows a more 



direct evaluation of the logarithmically-enhanced corrections to the heavy-quark total hadropro- 
duction cross-section, up to NLL accuracy^. A key simplification in this method is the direct 
diagonalization of the radiative factor A^jy f°r the resummed parton-level total cross-section in 
the colour-basis defined by colour- singlet and colour-octet final states. Details on the calculation 



and the illustration of the general method are presented elsewhere || [12[ . In this paper we limit 
ourselves to presenting the final theoretical results on soft-gluon resummation to NLL order, and 
we mainly concentrate on the ensuing phenomenological predictions. 



A similar structure also holds for the production of dijets at fixed invariant mass 1 1 . 

5 We have checked that our resummation formulae are consistent with the NLL results on the heavy-quark 
cross-section at fixed invariant mass obtained in ref. [pi. 
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The rest of this paper is organized as follows. In Section |2] we describe our notation, and 
collect the fixed-order leading and next-to-leading formulae that will be used in the following. In 
Section |3] we give the formulae for the NLL resummed cross-sections, and in Sect. f| we present 
our phenomenological results. 



2 Notation and fixed-order calculations 



In this Section we recall the known theoretical results [[y| 14| on the calculation of the heavy- 



quark hadroproduction cross-section at LO and next-to-leading-order (NLO) in QCD perturbation 



theory. We shall follow closely the notation of ref. 13 



The total cross-section a is given by the following factorization formula 

a(p h , m 2 ) = ^2 f dxi dx 2 F^Xt, p 2 ) Fj(x 2 , p 2 ) (p; m 2 , a s (p 2 ),p' 



(5) 



where m is the mass of the heavy quark, S is the square of the centre-of-mass energy, i, j denote 
parton indices (i = q, q, g) and the dimensionless variables p and ph are 

4m 2 p h 



Ph = —q- , P = • (6) 

b X\ x 2 

The parton densities F t (x, p 2 ) and the partonic cross-sections <7y depend on the factorization scale 
p (which we always set equal to the renormalization scale) and on the factorization scheme. We 
use the MS factorization scheme. 

We also introduce the dimensionless functions by rescaling the partonic cross-sections as 
follows 

cr l:j (p]m 2 ,a s {p 2 ),p 2 ) = ^iJL- fijfaasQjt 2 )^ 2 /™ 2 ) . (7) 
The functions are perturbatively computable and given by the following expansion in a s 



fij(p; a s (p 2 ),p 2 /m 2 ) = f-j{p) +Aira s {p 2 ) 



4"(p)+7!>)ln^ 



E. 

m 1 



J2<(p 2 )f^(p;p 2 /m 2 



n=2 



The LO terms ffj in eq. (^) are explicitly given by 



7T T R C F 
6 N r 



71 



12 JV C 2 - 1 



PP (2 + P) 



(3p\SC F 



(4 + 4p-2p 2 )i l n i±|-4-4p 



+ C; 



where (3 = \/l — p and flf 1 (p) = for all the other parton channels. 



3p 



21 lni±4-4-5p 







AD 



(9) 



(10) 



The NLO terms fl/(p) in eq. @ are obtained by simply performing the convolution of f-j (p) 

with Altarelli-Parisi splitting functions. The remaining NLO contributions fij\p), first evaluated 
in ref. , are not known in analytic form but only numerically |TB|, [TJJ . An analytic parametriza- 
tion of the numerical results is available in ref. JOJ. The functions fqg (p) and fj$(p) are plotted 
in fig. D (solid lines). 
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2.1 Threshold behaviour 



We are mainly interested in the behaviour of QCD corrections near the threshold region, p — ► 1. 

AO] 

lij 



In this region, the LO functions fjj\p) vanish because of phase-space suppression: 



^) P 3iW*^ ' /S'W^i^r^-^j^-o, (ii) 

while the NLO functions do not. The behaviour of the latter is analytically known and given by 
the following expressions [Tj| [TJj 

f&\p) + 72 (P) In ^ = ^ 0) (P) { (C F - ±<7 A ) g + 2C F ln 2 (8/? 2 ) (12) 



2/3 

(8CV + CU) ln(8/3 2 ) - 2C F ln(4/3 2 ) In ^ + C 2 ( ^ J + 0(1 - p) 

m z \ m I 



StM + 7>) m £ = ^/£>M { ^_ 2) $ + ^ '" 2 (83 2 ) (13) 

and fij\p), fij\p) = for all the other parton channels. 



The first term in the curly brackets of eqs. ( |l2"l , |T3l) is due to the Coulomb interaction between 
the heavy quark and antiquark. The remaining logarithmic contributions are produced by soft- 
gluon bremsstrahlung. Note that the right-hand sides of eqs. (|l|,|13|) also include the constant 
coefficients C 2 , C3, which are due to large-momentum virtual corrections. The explicit expressions 
for these coefficients are 



where riif denotes the number of massless flavours (i.e. the number of quark flavours that are 
lighter than the produced heavy-quark pair) and 6*2,3 are obtained from the values of the constants 
a reported in Table 1 of ref . El : 



768 

C 2 = 36vraf ~ 20.45 , C 3 = — ira g 9 ~ 37.23 . (15) 



In the strict threshold limit the Coulomb contributions dominate the radiative corrections in 
eqs. ( |T2"| , |T3"| ) . However, in most of the hadroproduction processes of phenomenological interest, the 
heavy-quark pair is not produced exactly at threshold. In these cases, it turns out [|13|, [TJ] that the 
logarithmic terms due to soft-gluon emission give the bulk of the NLO corrections. Multiple-gluon 
radiation at higher perturbative orders leads to more enhanced logarithmic corrections and the 
coefficient functions f^\p;p 2 /m 2 ) in eq. (§) behave as 

/IrWV^WfV)^ 2 • (16) 

Resummation of these soft-gluon effects to all orders in perturbation theory can be important to 
improve the reliability of the QCD predictions. 
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2.2 N- moment space 

The resummation program of soft-gluon contributions has been carried out JI], |2| in the Mellin- 
transform space, or TV-space. Working in TV-space, one can disentangle the soft-gluon effects in the 
parton densities from those in the partonic cross-section and one can straightforwardly implement 
and factorize the kinematic constraints of energy and longitudinal-momentum conservation. 

The Mellin transform of the heavy-quark hadronic cross-section a(ph,m 2 ) is defined as follows 



a N (m 2 ) = I dp p N 1 a(p,m 2 ) 



(17) 



where TV" is the moment index. In TV-moment space eq. ([5]) takes a simple factorized expression 



o N [ra 



F i)N+1 (p 2 ) F jtN+1 (p 2 ) a ijtN (m 2 ,a s (p 2 ),p 

a 2 s (p 2 ^ 



l :J 



53 Fi,N+i(p 2 ) F jjN+1 (p^ 



l :J 



fij,N{.ct s {p ),fi /m ) 



(19) 



where we have used eq. (|8|) and we have introduced the TV-moments ^^(/z 2 ) of the parton densities, 



F itN {p 2 



doc X ^ X j j-L ^ . 



(20) 



and likewise for any other x-space function. 

The TV-moments of the LO functions in eqs. (||) and ([H]) are given by the following explicit 
expressions [l7j 



AO) 

Jqq,N 



AO) 
J 99,N 



7T2 T R C F T(N + 1) 



4 TV C r(7V + 7/2) 



(TV + 2) 



(21) 



7T2 



T, 



R 



T(N+1) 1 



x 



4 TV C 2 - 1 r(7V + 5/2) TV + 3 

TV 3 + 97V 2 + 207V + 14 TV 2 + 87V +11 

2Cf t—z — . . — — — C a 



(7V + l)(7V + 2) 



27V + 5 



(22) 



To discuss the structure of the higher-order corrections considered in the next Section, it is 
useful to decompose the various cross-section contributions according to the colour state (singlet 
1 and octet 8) of the heavy-quark pair produced in the final state. When applied to the LO 
functions /A- , this decomposition gives 



(o) 



with 



(23) 



f (0) (o) = ■ 

Jgg,xKP) Nc (]S[*-l) 8 



Tf, 



TTpp 



(4 + 4p - 2p 2 ) - In 



1 , 1 + 



-4-4p 



(3 1-/3 

4- } i(p) = o, &( P ) = fS\p) 



(24) 
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The iV-moments of /i?i(p) are 



3 

7T2 



AO) _ 

J 99,1, N — A 



T(N + l) A^ 3 + 9iV 2 + 2QN + 14 



4 iV c (iV c 2 -l) r(A^ + 5/2) (iV+l)(iV + 2)(iV + 3) ' 



(25) 



and those of ff)${p) an d s(p) are then obtained from eqs. by using eq. (^| 



e(0) 



We also introduce the functions fjj ° onl (p) that describe the Coulomb contributions to the 



NLO corrections in eq. ([H],|nj). Using the singlet/octet colour decomposition, we have 



iS Coul (P) = mp)c F 



(0) 



'■I.) 



i/9 



iS Coul (p) 



whose Mellin transforms are given by 

f (l),Coul 7T T R Cf 

Jqq,S,N ~ 



3N + 5 



,(l),Coul 

Jgg$,N 



,(l),Coul 
J 99,1, N 



96 A" c 2 (N + l)(N + 2) 

Tfl 



(26) 
(27) 

(28) 



+ 3 



192 N C {N% - 1) 
N 2 -2 



N r . ( 3I N+2 - 



64 N* 



21^ + 2In + i — In+2 



2/jv + 2/jv+i _ In+2 — 



N+l N+2 
2 2 



N+l N+2 
2 2 



N+l N+2 



The function I N is defined by: 



in 



We can rewrite In as follows 



in 



dp p 



N 



dp p 



1 



N 



1 + VT^P 



dn 



In - + 21n(l + Jl-p) 
P v 



(29) 
(30) 

(31) 

(32) 



The integral of the first term can be done exactly, while to integrate the second term we first 
expand the logarithm of 1 + y/1 — p as: 

4 

ln(l + x) = ^2 a n x n 



(33) 



71=1 



with a n = 0.9991, -0.4828, 0.2477, -0.0712 for n = 1, ... ,4, and then perform the integral term 
by term. The expansion in eq. ([33]) is accurate to 10~ 4 in the needed range < x < 1. The final 
approximate result for I N is: 



IN 



T(N + 1) T(N +l)T(^) 

* , ,\ MN + 3/2) - tp(N + l)] + 2Ya n ^— — 2 A 

T(iV + 3/2)^ v ' ' V n ^ T(N + n/2 + 3/2) 



(34) 



The threshold region p — > 1 corresponds to the limit N — > oo in A^-moment space. In this 
limit the Mellin transform of the threshold expansions C[l2| , p~3f) of the NLO corrections can be 
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easily computed. Having introduced the Coulomb contributions /i 1 ° ul , we can write the result 
as follows 



/?g,iV + Jqq,N m m 2 



1 f (0) 
Jqq,N 



An 2 



2C F \n 2 N+ [C A + AC FlE + 2C F In 



4m 2 



In AT 



+ C qq ^/m 2 ) + 0^ + f ( 



(0) 



2C A In 2 AT + 



N 2 



N 2 



(l),Coul 
qq,8,N ' 

!c A + 4C A7 £ + 2C A ln /; 



(35) 



4m 2 



In AT 



(l),Coul 
99,8,^ ' 



(36) 



where ^ E = 0.5772 ... is the Euler number and the constant coefficients C qq , C gg are related to 
those in eqs. (|l4l,|T5|) : 



Cqq{ 


= C 2 (^- 2 ) + 2C F 
m z 

+ (8C F + C A )( lE 


In 2 2 + ( 7J5 - 
-2-ln2) , 


) 4m 2 


vr 2 

+ y 


+7^(7^ 


-4) 


Cgg ( o ) 

m 2 


= C 3 (^- 2 ) + 2C A 
m z 


ln 2 2 + ( 7i? - 


2 ) ln f% 

Am 2 


7T 2 

+ y 


+ 7s (7s 


-4) 




+ c ^" 2 20) (-- 2 -'-) 











(37) 



(38) 



The double and single logarithmic terms In 2 A" and ln A" in eqs. ( |3~5"1 , |3~6"D are produced by soft- 



The 



gluon radiation. These terms dominate the corrections in the curly brackets of eqs. (p5|,p 
resummation of soft-gluon contributions to all orders in perturbation theory is considered in the 
following Section. 



3 Soft-gluon resummation 

The partonic cross- sections ct^jv with ij 7^ qq, gg vanish at LO and are hence suppressed by 
a factor of a s with respect to & q q t N and o" 99 ,at. In the threshold or large- A" limit, this relative 
suppression is further enhanced by a factor of 0(1/N). Therefore, we make no attempt to resum 
soft-gluon corrections to these partonic channels and we restrict ourselves to & q q t N and <J gg ,N- 



3.1 Total heavy-quark cross-section to NLL accuracy: the radiative 
factors 

As discussed in the Introduction, the resummed radiative factor [A^jvJu in eq. ([|) depends 
in general on the partonic subprocess (labelled by ij) and on its colour states (labelled by I, J). 
However, in the case of the threshold behaviour of the total cross section, most of the complications 
related to the colour structure can be easily overcome. Up to NLL accuracy, it turns out |], |12| 
that the colour matrix A^jy is diagonal with respect to the basis I = {1, 8}, 1 and 8 being the 
singlet and octet state of the produced heavy-quark pair. We can thus present the exponentiated 
formulae for its eigenvalues Aj^jy in the two colour channels. 
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We use a notation similar to that in ref. to facilitate the comparison with the known results 
for the DY process. In the MS factorization and renormalization schemes, the NLL resummed 
expressions for the radiative factors A%j,i,N are: 



In A 



qq,l,N 







In A 



qq,8,N 



i z N ~ 1 -l 

dz 

1 — z 

N-l _ i 



2A q (a s {ql))+0{a s {a s \nNY 



dz 



z 



1 



z 



Ql 



2A q (a s (ql)) + D QQ (a s (4m 2 



'1-z) 



+ 0{a s {a s \nNf 



In A 



99A, N 



In A 



gg,s,N 



dz 



dz 



z 



N-l 



- 1 



1-Z 
N-l _ ^ 



4nv»(l-*)* rf g 2 

■? qj_ 



0{a s {a s \nN) } 



z 



z 



W (1 _ 2) 2 dg 2 



2A g (a s (q] 



D 



QQ 



(a. (4m 2 (1 - z)< 



+ 0(a s (a s \nN) k ) 



where A q ]15j, A g [pT6| , -Dqq [§] are the following perturbative functions 



A q 




— Op — 

7T 


1 + — - 

. 2 7T_ 


+ , 


Ag 


>*) 


= C.Hl 

7T 


1 + 1 

2 7T . 


+ 0{al) , 


Dqq 




= — — 

7T 







and the coefficient K is given by: 

K = 

Some comments are in order. 

The LL terms a™ ln n+1 N in eqs. 
function Dqq and by truncating 
resummed results first derived in ref. 



67 
18 
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(39) 

(40) 
(41) 

I 

(42) 

(43) 

(44) 
(45) 

(46) 



2D are obtained by neglecting the contribution of the 
sj and A g (a s ) to their first order. One thus recovers the LL 



The contributions of the functions A q ,A g , Dqq have a direct physical interpretation. The 
function A q (A g ) measures soft and collinear radiation from the qq (gg) partonic channel in the 
initial state. This is the only contribution that appears in hadroproduction processes of colourless 
heavy-mass systems, such as the DY process [fT], p|, |TT|] or Higgs boson production via gluon 
fusion |L6|, [18 . 



The function Dqq is due to soft emission from the final-state heavy-quark pair. Since the 
final-state quarks are massive, they do not lead to collinear logarithms and the function Dqq only 
enters at NLL accuracy. At threshold, the heavy-quark pair acts as a single final-state system 
with invariant-mass squared Q 2 = 4m 2 and with colour charge given by the total QQ charge. 
Thus, its contribution to the radiative factors A^jv vanishes in the colour-singlet channels (see 
eqs. ( |3l3|j4"l"D ) and is proportional to Ca and independent of the initial state in the colour-octet 
channels (see eqs. ( f4T|f42D ) . 

Since we know the resummed radiative factors Ay i # only to NLL accuracy, we use the 
expressions (|39"l - f42"D by explicitly carrying out the z and q]_ integrals and neglecting terms beyond 
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NLL. Thus, we write: 

iN (a a (ji 2 ),fi 2 /m 2 ) = exp{hiN g§\boa a (n 2 )hiN) + gf^%a s {^ 2 )\nN^ 2 /m 2 ) 

+ 0(a 2 s (a s In N) k )\ . (47) 



The functions and gf( 2 ) resum the LL and NLL terms, respectively. They can be explicitly 
computed as in refs. p], 0. We find: 

= ^TT PA + (1 - 2A) ln(l - 2A)] , (A) = <$(A) , (48) 



iqqV> tt6 A l ' v 1 K VJ ' »" v> C F * qq 



and 



4i A. H m) = -7s— 7- In 1 - 2A + — - - 

7TOo VTOg 



2A + ln(l - 2A) + - ln 2 (l - 2A) 



' A [2A + ln(l - 2A)] - ln(l - 2A) In ■ P 



2tt 2 6 2 ) n vrfeo 4m 2 ' 

»Sx(A,Ai 2 M a ) = ^£!i(A,/i 2 /™ 2 ), (49) 



^ 2) 8 (A,/x 2 /m 2 ) = ^g) 1 (A, /U 2 /m 2 )-^-ln(l-2A). 

where &o> &i are the first two coefficients of the QCD /^-function 

= \\C A -4T R N S = Y1C\ - 10C A T R N f - QC F T R N f 

12tt ' 1 24tt 2 ' [ } 

Note that the LL function g^ N depends on the initial-state partonic channel, but it is independent 
of the colour state. 

It is common practice to study the scale dependence of the cross-section in order to assess the 
impact of neglected higher-order effects. Normally, in the fixed-order calculations, one considers 
the range m/2 < fi < 2m. When the threshold region is approached, however, a new scale comes 
into play, namely the typical distance from threshold of the hard production process. Thus, the 
process is really a two scale process, and one is forced to use a wider scale range in order to estimate 
the error. When the NLO calculation is improved by LL resummation, this is no longer needed, 
since the large threshold logarithms have been properly resummed. Thus, one can still safely 
consider the range m/2 < fi < 2m. Observe, however, that the inclusion of LL resummation does 
not bring a reduction of the scale dependence in the range m/2 < fi < 2m. The improvement comes 
simply from the fact that in this case we are allowed to use scales of order m. On the other hand, 
when NLL resummation is used, one does expect a reduction of the scale dependence in the range 
m/2 < n < 2m. This can be seen explicitly from eqs. (f4T[ - f49|) . A scale variation from fi to \i! is 
perturbatively equivalent to the replacement a s (/z 2 ) — > a s (fi' 2 ) = a s {fi 2 )[l — b a s ([i) ln// 2 /yU 2 + . . .]. 
This replacement in the LL functions gQ (boa s In N) leads to additional NLL terms in the exponent 
on the right-hand side of eq. (f|7|). These additional terms are partly cancelled by those that 
originate from the explicit /i-dependence of the functions glj\(boa s In N, fi 2 /m 2 ). The left-over 
NLL contributions exactly match the scale dependence of the parton densities at large values of 
N. For instance, in the qq channel we have 



A OT -i,7v la s (/i /2 ),^j =A q g >IiN U s (/i 2 ),£ I jexp|2^a s (/i 2 )lnA^ln^ + C'(a s (a s lnAr)' 



(51) 
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and, in the calculation of the full hadronic cross-section, the exponential on the right-hand side is 
compensated by the scale dependence of the quark densities in the MS factorization scheme: 

F q , N {fj! 2 )F q>N {fj,' 2 ) =F^(/i 2 )F^( / i 2 )exp|-2^« s (^ 2 )lniVln^ + 0(« s (« s lniV) /£ ) j . (52) 



3.2 NLL resummed cross-section 



Using the results of Sect. [3J], we can introduce the NLL resummed partonic cross-sections as 
follows 

4 e £V,«^ 2 )V) = ^ E fS:A(^ 2 ),^ 2 /m 2 ) , (53) 



1=1,8 



f^i%M^ 2 ),^ 2 /m 2 ) = 4-i^(a s (/i 2 ),/i 2 /m 2 ) A ij>ItN+1 f « s (/i 2 ), £^ J • 



(54) 



Note the mismatch between the moment indices in the right-hand side of eq. fl5"4"|). The radiative 
factors Ay^jv+i depend on the moment index N+l, like the parton densities in the factorization 
formula (|T5|). 



The terms fi^N m ec L- 



are given by the the LO functions in eqs 



after 

correction by the Coulomb contributions in eqs. (|28|j29|j30| ) and by the iV-independent coefficients 
Cy in eqs. (0,0) 



f^ = {fSl N +^ s ^)fSi 



Coul 
N 



7T 



-C^/m 2 ) 



(55) 



Using the definition in eq. ( |5"5] ) and the perturbative expansions of the radiative factors A^ati 



l n ^qq,l,N 
In ^ij,8,N 



9* 



In A 



ggA,N 



71 



2C F In 2 N + {AC fIe + 2C F In^/m 2 ) IniV + 0{a 2 ) 



\nA ijiliN + ^^C A lnN + 0{a 



2 ) 

s) i 



(56) 



it is straightforward to check that eqs. (0,0) correctly reproduce the NLO threshold behaviour 
in eqs. (0,0). In particular, since both the singlet and octet states are produced at LO via 
gg fusion, the colour factor {N 2 — A)/{N 2 — 2) that appears in eq. (0) is due to the fact that 
final-state soft-gluon radiation is weighted by the ratio 



(0) 

AO) 



N 2 



N 2 

c 



[1 + (D(1/N)) , 



(57) 



which differs from unity near threshold. 

The all-order factorization of fijY^ with respect to the radiative factors Ajj jIjA r +1 in the re- 
summed partonic cross-sections (Q) is justified by the fact that the 0(a s )-terms in eq. (0) are 
produced by non-soft virtual corrections to the LO subprocesses. The resummation of Coulomb 
corrections was considered in ref. and it turned out that, in the experimental configurations 
of practical interest, its quantitative effect is much smaller than that due to soft-gluon radiation. 
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Because of this reason, in our resummed formulae we do not include pure Coulomb effects beyond 
NLO, but simply the soft-gluon corrections to the 0(a^) Coulomb contribution. 

Note, however, that we include the constant terms CV,(/i 2 /m 2 ) as in eqs. ( |5^J55| ). This proce- 
dure is analogous to that used in NLL resummed predictions for other observables, such as e + e~ 
event shapes []r9 |. As shown in Sect. O] , these constant terms are important to accurately match 
the exact near-threshold behavior at NLO. Moreover, when combined with the NLL soft-gluon 
factors Aij^ N , they lead to the following expansion 



l + ^-C^/m 2 ) 



TT 



A ijXN (a s (/j 2 ),/j 2 /m 2 ) = 1 + ofib 1 



hi 



2n 



N + c nM . x In 2 "" 1 N 
- c„ i2n _ 2 (/i 2 /m 2 )ln 2 "- 2 N + 0(ln 2n ' 3 N)] ,(58) 



n=l 



that correctly predicts the value of the coefficients c„ j2 „_ 2 (/x 2 /m 2 ), as well as of c nj2 „ and c n>2n -i 
that are independent of Cy. Note also that, while c„ )2n and c„ >2 „_i are \i- independent, the coeffi- 
cients c n>2 „_ 2 (/i 2 /m 2 ) are not. Their dependence on /i is obtained by combining that of CV,(/z 2 /m 2 ) 
with the explicit scale dependence of A i jj ) Ar(a s (/i 2 ), /i 2 /m 2 ) at NLL order. 

Including the Cy constant term can also be viewed as using the known NLO cross-section to 
improve the resummation formula beyond the NLL approximation. This can be seen schematically 
in the following way. We write the left-hand side of eq. fl58" ) as: 



In [(1 + a s /nC)A] = In A + a s /nC , 
and since our NLL expression for In A has the form 

In A = In Ng {1) (a s In N) + g {2) (a s In N) 



(59) 



(60) 



according to eq. (§), the term a s /nC corresponds precisely to the 0(a s ) term in the expansion of 
the NNLL contribution a s g^(a s In N). 

We use soft-gluon resummation to NLL accuracy to introduce an improved hadronic cross- 
section follows 



a 



(res) 
N 



m 



a 



(res) 
ij,N ' 



m 2 ,a s (/x 2 ),/i 2 ) - (&§$(m 2 ,a a (p 2 ), 



a 



(NLO) 
N 1 



m 2 ) 



where a\^ L °^ is the hadronic cross-section at NLO, (rj^y is given in eq. 



and (<7 



. (res) 

hj,N 



(61) 



rep- 



resents its perturbative truncation at order aP s . Thus, because of the subtraction in the square 
bracket on the right-hand side, eq. ([H]) exactly reproduces the NLO results and resums soft-gluon 
effects beyond 0{a\) to NLL accuracy. This defines our NLO+NLL predictions. 

The resummed formulae presented so far are given in iV-moment space. To obtain cross-sections 
in the physical x-space, one has to perform the inverse Mellin transformation: 



a 



(res) 



(p,m 2 ) = -i T / dNp- N a^' 



m 



(62) 



When the iV-moments <jn are evaluated at a fixed perturbative order in a s , they are analytic 
functions in a right half-plane of the complex variable N. In this case, the constant Cup that 
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defines the integration contour in eq. (|62]) has to be chosen in this half-plane, i.e. on the right of 
all the possible singularities of the iV-moments. 

An additional complication occurs when the iV-moments are computed in resummed pertur- 

(*), 



bation theory. In this case, since the resummed functions ghi(A) in eqs. (H^,§^) are singular 



at A = 1/2, the soft-gluon factors A i jj tN (a s (jj 2 ), jj 2 /m 2 ) in eq. QI7D have cut singularities that 
start at the branch-point N = N L = exp(l/26 a s ). These singularities, which are produced in 
eqs. ( |39l - [42] ) by the g_i_-integration down to the Landau pole of the running coupling a s (q 2 _), signal 
the onset of non-perturbative phenomena at very large values of iV or, equivalently, in the region 
very close to threshold. 

The issue of how to deal with the Landau singularity in soft-gluon resummation formulae 
for hadronic collisions was discussed in detail in ref. 0]. In the evaluation of the inverse Mellin 
transformation (02) we thus use the Minimal Prescription introduced in ref. The constant 

Cup is chosen in such a way that all singularities in the integrand are to the left of the integra- 
tion contour, except for the Landau singularity at N = Nl, which should lie to the far right. 
This prescription is consistent |7]] with the perturbative content of the soft-gluon resummation 
formulae because it converges asymptotically to the perturbative series and it does not introduce 
(unjustified) power corrections of non-perturbative origin. These corrections are certainly present 
in physical cross-sections, but their effect is not expected to be sizeable as long as m is sufficiently 
perturbative and p is sufficiently far from the hadronic threshold. Obviously, approaching the 
essentially non-perturbative regime m — > 1 GeV, p — > 1, a physically motivated treatment of non- 
perturbative effects has to be introduced. In the following Section, we limit ourselves to present 
numerical and phenomenological results for kinematic configurations in which non-perturbative 
corrections should be smaller than the estimated uncertainty of the perturbative predictions. 

In all cases considered, we have verified that the minimal prescription is in practice equivalent 
to evaluate eq. (|62|) order by order in perturbation theory, until the corrections become numerically 



insignificant. An example of this procedure will also be shown in Section [O 



4 Results 

We present in this Section some numerical results, to provide an illustration of the size of the 
effects considered and a new estimate of heavy-quark production cross-sections of relevance for 
present and future hadronic colliders. 



4.1 Parton-level results 

We start by discussing the resummation effects at the level of partonic cross sections. The re- 
summed partonic cross-section for the production of a heavy quark pair can be obtained from 
eqs. ( pT]) and (p2"| ) by assuming structure functions of the form F(x) = 5(1 — x), and therefore 
Fn = 1 for all complex values of N. 

We consider first the 0(af) terms in the expansion of the resummed cross-section, in order 
to estimate to which accuracy this reproduces the exact NLO results. In fig. [1] (left) we plot the 
function fj^\ defined in eq. (H), as a function of r] = (1 — p)/p. The exact 0(a^) result [13| 



is 



compared with three possible implementations of the resummation procedure, all equivalent at 
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Figure 1: Left (Right): the function fqq(p) (fgg\p)), plotted as a function of rj = 
(1 — p)j p. The solid line represents the exact NLO result |l3|, the short-dashed line 
corresponds to the O(a^) truncation of the resummed result defined by eqs. (]53|-[55|); the 
dot-dashed line is obtained from this last result by setting the constant C qq (C gg ) to 0; 
the dashed line is obtained instead by the replacement in eq. (]63|), with A = 2. 



NLL. In one case (dot-dashed line) we set the constant C qq introduced in eqs.fl3~7|) and (^) equal 

to 0. In the second case (short-dashed line) we include the contribution of C qq . In a third case 

we correct the contribution of the constant C qq by a term which is suppressed by a factor oll/N, 

which does not introduce poles on the real N axis, and which gives vanishing first moment to 
f (i) . 

Jqq,N- 

Cij -> Qj (1 - N _ 1 ) , ij = qq, gg . (63) 

In our applications we shall consider the two cases with A = (namely no correction to the 
term) and A = 2 as a way to establish the size of subleading corrections beyond the NLL order. 
In all cases we include the effect of the leading-order Coulomb effects, as described in eq. (|55|). 

As one can see from fig. [I], the inclusion of the finite term C qq is essential to accurately 
reproduce shape and normalization of the function f q ^ near threshold. The agreement deteriorates 
unavoidably for /> > 1, as here terms subleading in 1/N become important. Note that the 
two choices A = and A = 2 braket the exact result, and thus provide a good estimator of 
the subleading terms' systematics. The choice A = 2, furthermore, provides a very accurate 
description up to values of r] of order 1. This is the region which dominates the production 
cross-section in the cases of interest. 

Analogous results for the gg channel are given in the right panel of fig. |I[ The agreement is 
again very good near threshold. Far above threshold the exact NLO result is dominated by the 



t-channel gluon exchange diagrams [0, 20 1, which give rise to a 1/N pole not controlled by the 
soft-gluon resummation. 

The fully resummed parton-level cross-sections are shown in fig. ^| for the qq and gg channels 
(left and right panel, respectively). Here and in the following we shall define the resummed 
cross-sections as in eq. (|BTD, that is, we substitute their O(a^) terms with the exact NLO result, 
using the same choice of renormalization scale /i. In this way our results are exact up to (and 
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v=i l -p)/p v=( 1 ~p)/p 



Figure 2: Partonic cross-section for the processes qq — ► QQ (left) and gg — > (right) 
(in pb, and for mg = 175 GeV). The dashed line is the exact NLO result [L4|]; the 
short-dashed (solid) lines correspond to the NLO+NLL result, with the coefficient A 
defined in eq. ( p3| ) equal to (2). The lower and upper curves correspond to inclusion or 
neglect of the Coulomb contribution. 



including) O(a^), and include the NLL resummation of terms of O(otf) and higher. We compare 
the fixed-order results (dashed lines) with the resummed results. For these we provide both the 
A = and A = 2 prescriptions, as well as the inclusion and neglect of the Coulomb terms in 
eq. fl55|) . Note that even at the level of resummed cross-sections the prescriptions A = and 
A = 2 braket the fixed order result. Note furthermore that, for the curves without the Coulomb 
contribution, the oscillations of the cross-section near threshold are significantly reduced relative 
to the LL resummation which was presented in fig. 3 of ref. |7|]. 



4.2 Hadron-level results 

In this Section we present results for the full hadronic cross-sections. As a default set of parton 
densities, we shall use the MRSR2 set described in [2lj|. For the top-quark mass we shall use, 
unless otherwise indicated, m t = 175 GeV. 

Figures || f| and |5] present the ratios: 

^-res _ .-.NLO ~res _ „NLO ^res _ -.NLO 

NLL,gg °gg °NLL,qq °qq ^NLL.feg+qq) ° (gg+qq) 

^NLO ' ^-NLO ' rr NLO ' ^ ' 

u ss u qq U (gg+qq) 

For each channel we present the results using both the A = (dashed lines) and A = 2 (solid lines) 
prescriptions. We also show the dependence on the choice of renormalization and factorization 
scales, which we always take equal, and varying within the set /i = (m t /2, m t , 2m t ). Note that the 
size of the resummation effects is larger for the larger scales, contrary to the behaviour of the scale 
dependence of the NLO cross section. This suggests that the scale dependence of the resummed 
cross-section will be reduced relative to that of the NLO results. The size of the resummation 
corrections is large at small centre-of-mass energies, as should be expected for the region where 
production is dominated by the threshold region. At Tevatron energies (VS= 1.8 TeV), the size 
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Figure 3: Contribution of gluon resummation at order 0(aj.) and higher, relative to 
the exact NLO result, for top-pair production via qq annihilation in pp collisions. The 
solid (dashed) lines correspond to A = 2 (A = 0). The three sets of curves correspond to 
the choice of scale fi = 2mt, nrit and mt/2, in descending order, with rrit = 175 GeV, and 
PDF set MRSR2. 
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Figure 4: Same as fig. 0, for production via gg annihilation. 
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Figure 5: Same as fig. ||, for the combined production channels gg + qq. 

of the corrections equals 4% (6%) for A = 2 (A = 0) and ji = m t . It becomes equal to 0.4% 
(0.5%) for /i = mt/2 and to 9% (14%) for fx = 2m t . 

The scale dependence of the resummed cross-section, compared to the NLO one, is given in 
figs. ^| and 0, corresponding to the A = and A = 2 cases, respectively. Note the significant 
reduction in scale dependence, more marked in the A = case. More importantly, note that the 
band of variation of the resummed cross section lies entirely within the band of variation of the 
NLO cross-section for Tevatron energies and above. This shows that previous estimates of the 
theoretical systematic uncertainty for the Tevatron cross-section were correct, and can now be 
improved thanks to the NLL calculation presented here. 

To display the importance of the inclusion of the NNLL CV,- terms, we show the same plot of 
the scale dependence with CV, = in fig. |8|. While the scale sensitivity is slightly worse than in the 
cases with CV, 7^ 0, there is still an important improvement over the NLO behaviour. In the same 
figure, we also show the effect of neglecting the contribution of Coulomb terms of order O(af) and 
higher. With the exception of the low-energy points, where Coulomb effects are very important 
because of the closeness of the threshold, inclusion of the Coulomb terms in the resummation is 
not a significant effect. 

A more detailed representation of the scale sensitivity of the resummed top cross-section is 
shown in fig. |9|, where we show the scale dependence in the wide range m t / 10 < ft < lOm^. 

It was pointed out in ref. that large NLL contributions to the resummation of leading 
logarithms arise from the inclusion of NLL corrections to the relation 1 — z N ~ 8(1 — z — 1/-/V), 
which can be used to perform the integrals in eqs. (|39D-([4~2j). This relation, which is valid at LL 
level, should be replaced at NLL order by 1 — z ^ 9(1 — z — e~ lE /N). Neglecting terms beyond 
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Figure 6: Scale-dependence of the ti production cross-section in pp collisions, as a 
function of y/S. The solid lines represent the exact NLO result for different choices of 
\x {ji = mt/2 and 2mt), normalised to the \i = mt result. The solid lines represent the 
NLO+NLL result (with ^4 = 0) for different choices of jx ([A = mt/2, mt and 2mt), 
normalised to the NLO \i = mt result. 
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Figure 7: Same as fig. |, but for A = 2. 
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Figure 8: Same as fig. || but with the NNLL terms proportional to set to zero 
(short-dashed lines). The dashed lines represent the same calculation, with the Coulomb 
contributions in eq. (|55|) removed. 
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Figure 9: Scale dependence of the total tt production cross-section in pp collisions at 
y/S = 1.8 TeV. LO (short-dashed line), NLO (dashed line) and NLO+NLL (solid and 
dot-dashed lines, for A = 2 and A = respectively). 
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Figure 10: Same as fig. |6|, but with the NLL result replaced with the resummation of 
the leading logarithms only, and inclusion of the terms defined by eq. (|65|). 



NLL accuracy, this replacement amounts to the substitution^: 

A "^( 1 + ^ai^v) A "- < 65 > 

For the specific choice fx = m t , we showed in [[jj that the 0(a^) truncation of the LL resummed 
calculation provides a better agreement with the exact NLO result for the partonic cross-section 
when this substitution is applied. This agreement is however accidental and limited to this partic- 
ular choice of the renormalization scale. Figure [10] shows in fact that, with the sole substitution 
in eq. (|65"D , no improvement in the scale dependence is observed relative to the NLO calculation 
when the scale is varied. Such an improvement can only be obtained when the full set of NLL 
terms is included. 

The heavy-quark pair cross-section in pp collisions at y/S = 1.8 and 2 TeV are shown, as a 



function of the heavy-quark mass, in figs. [11] and [12]. We chose the value A = 2, which gives 
a more conservative estimate of the scale uncertainty, and PDF set MRSR2 (a s (Mz) = 0.119). 
The two dashed lines correspond to the NLO result, with fx = m t /2 and 2m t . The two solid lines 
correspond to the fully resummed NLL result. At \/S = 1.8 TeV and for m t = 175 GeV, we get a 
fully resummed result of a t t = 5.06^°'^ pb, compared to the fixed NLO one of a t i = 4.87^°'^ pb 
(the central values representing the results with fx = m t ). The scale uncertainty of the resummed 
result is reduced by almost a factor of 2. The current experimental results from CDF |22|] and 



DO |3| are respectively: <r a = 7.6 +j ° pb (CDF, at m t = 175 GeV), and a a = 5.5 ± 1.8 pb (DO, 



at m t = 173.3 GeV). 



6 This is also equivalent to adding to lnAfl i the sole contributions proportional to 'je in the expressions (k 



for 5l (2) . 
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Figure 11: Total tt production cross-section in pp collisions at y/S = 1.8 TeV, as a 
function of the top-quark mass. Dashed lines: NLO result; solid lines: NLO+NLL result. 
Upper lines: \i = mt/2; lower lines: fi = 2mt- 



The NLO and fully resummed cross-sections for LHC (pp collisions at \/S = 14 TeV) are 



shown as a function of the heavy-quark mass in fig. 113 



In Table [3] we collect the values of the cross-sections for m t = 175 GeV at the Tevatron and 
LHC. The rapid convergence of the higher-order corrections is displayed in table EL for the case of 
top production at the Tevatron. The sum of all entries in each row corresponds to the minimal 
prescription, while each fixed order term does not have any ambiguity due to the choice of the 
contour for the Mellin transformation in eq. (|62|). This supports the validity of the minimal 
prescription, since the truncated resummed expansion converges to it very rapidly. 

Potentially large resummation effects should be expected in bb production at fixed target 
energies, due to the closeness of the threshold. As an example, we collect in Table |] the values 
of the bb total cross-section calculated using the MRSR1 set (a s (Mz) = 0.112) for pp collisions at 
v^S* = 39.2 GeV, the configuration of the upcoming HERAB experiment at DESY. We show both 
the NLO results and the fully resummed ones, for several values of the bottom quark mass and for 
the three standard choices of renormalization scale. The resummation corrections are large and 
positive for large values of the renormalization scale, but become small - and eventually negative - 
for the smaller values of the renormalization scale. As a result, the overall scale dependence of the 
cross-section is significantly reduced in the fully resummed calculation, and the overall uncertainty 
band is fully included in the uncertainty band of the NLO calculation. 

The NLO cross-section, and the resummation effects, are largely enhanced by the choice of 
the MRSR2 set, for which the value of the coupling constant is larger. The results are shown in 
Table f|. The conclusions regarding the improved stability of the resummed result, as well as the 
reliability of the uncertainty estimate obtained from the scale variation of the NLO calculation, 
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Figure 12: Total ti production cross-section in pp collisions at = 2.0 TeV, as a 
function of the top-quark mass. Dashed lines: NLO result; solid lines: NLO+NLL result. 
Upper lines: fi = m t /2; lower lines: (i = 2m t . 



A t R = A t F 


pp at ^/S = 1.8 TeV 


pp at y/S = 2 TeV 


pp at ^/S = 14 TeV 


NLO NLO+NLL 


NLO NLO+NLL 


NLO NLO+NLL 


m t /2 


5.17 5.19 


7.10 7.12 


893 885 


m t 


4.87 5.06 


6.70 6.97 


803 833 


2m t 


4.31 4.70 


5.96 6.50 


714 794 



Table 1: Total tt cross-sections (m t = 175 GeV) at the Tevatron and LHC, in pb. PDF 
set MRSR2. 
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Figure 13: Total QQ production cross-section in pp collisions at y/S = 14 TeV, as a 
function of the heavy-quark mass. Dashed lines: NLO result; solid lines: NLO+NLL 
result. Upper lines: /j, = rriQ/2; lower lines: \i = 2rriQ. 





a 2 s 


a 3 s 


a 4 s 


a 5 s 


at 


«F 


qq 


3590 


766 


60.2 


2.4 


-0.6 


-0.08 


99 


298 


264 


98.1 


26.1 


5.8 


1.5 


qq + gg 


3888 


1030 


158 


28.5 


5.1 


1.4 



Table 2: Contributions to the total tt cross-sections (in fb) at the Tevatron (1.8 TeV) 
from higher orders in the expansion of the NLL resummed result, with /j, = mt, mt = 175 
GeV and PDF set MRSR2. The second column gives the exact NLO result. The last row 
only includes the sum of the qq + gg channels. 
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m 6 =4.5 GeV 


m 6 =4.75 GeV 


m fe =5.0 GeV 


NLO NLO+NLL 


NLO NLO+NLL 


NLO NLO+NLL 


m b /2 


27.4 27.3 


17.6 17.7 


11.2 11.5 


m b 


15.7 22.8 


9.90 14.9 


6.26 9.75 


2m b 


8.74 18.0 


5.47 11.7 


3.43 7.63 



Table 3: Total bb cross-sections (in nb) at HERAB (pp at = 39.2 GeV), as a function 
of the the b mass m b . PDF set MRSR1. 



/iR=/i F 


m 6 =4.5 GeV 


m 6 =4.75 GeV 


m fe =5.0 GeV 


NLO NLO+NLL 


NLO NLO+NLL 


NLO NLO+NLL 


m b /2 


45.2 44.9 


28.6 28.9 


18.0 18.6 


m b 


22.2 35.3 


13.9 22.9 


8.71 14.9 


2m b 


11.2 26.1 


6.94 17.0 


4.31 11.0 



Table 4: Total bb cross-sections (in nb) at HERAB (pp at y/S = 39.2 GeV), as a function 
of the the b mass m b . PDF set MRSR2. 
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Figure 14: Scale dependence of the total bb production cross-section at HERAB, at LO 
(dotted line), NLO (dashed line) and NLO+NLL (solid line). 



are nevertheless unchanged. A plot of the scale dependence of the LO, NLO and resummed 
cross-sections, for wif, = 4.75 GeV and PDF set MRSR2 is shown in fig. [14. 

As a comparison, we give here the value of the bottom-pair cross-section at the Tevatron, 
where we expect much smaller resummation effects. The NLO cross-section at the Tevatron (m& = 
4.75 GeV, MRSR2) is: o = 74.8, 56.2, 46.6 fib for \i = to&/2, m&, Im^. Including resummation we 
obtain: a = 74.2, 57.6, 52.2 fib for fi = rrib/2, mj, Im^. Some reduction in the scale dependence 
can be seen, but not at the level exhibited by the HERAB results. This is consistent with the 
fact that b production at the Tevatron takes place far away form the threshold, and the class of 
logarithms resummed in this work cannot be expected to give a dominant contribution. 



To conclude, we show in fig. [15] the effect of resummation on the total cross-section for fixed- 
target hadropro duct ion of charm-quark pairs. The typical values of 4m 2 / S at fixed-target energies 
(where \fS is in the range 20-40 GeV) are close to the values for top-pair production at the 
Tevatron. Contrary to the top case, however, charm production is dominated in pp collisions by 
the gg annihilation channel, due to the smallness of the antiquark sea. Considering in addition 
the larger size of the strong coupling constant at scales of the order of the charm mass, we 
therefore expect the resummation corrections to be large. Results are plotted as a function of 
the beam energy, for proton-induced reactions. The interval for the scale variation is limited to 
the range m c < fi < 2m c , since the parton densities we use (MRSR2) are limited to the domain 
fi 2 > 1.25 GeV 2 . Once more, the resummed result shows a significant reduction in the scale 
dependence over the NLO result. 
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Figure 15: Total cc cross-section in fixed-target pp collisions, as a function of beam 
energy, at the NLO and NLO+NLL. 

5 Discussion and conclusions 

We presented in this paper the first calculation of the heavy-quark hadro-production total cross- 
section accurate up to, and including, next-to-leading threshold-enhanced logarithms, resummed 
at all orders of perturbation theory. One can show that our analytical results for the resummed 
total cross-section are consistent with those obtained in in |J for the cross-section at fixed invariant 
mass of the heavy-quark pair. When the resummed partonic cross-sections are convoluted with 
hadronic parton densities, a significant improvement is observed in the stability of the results 
with respect to changes of the renormalization scale. In the case of top-quark production at the 
Tevatron, the corrections relative to the NLO calculation are of the order of 9% for \x = 2mt, 
4% for \x = rrit, and 0.5% for fi = m t /2. The result for \x = rrit can be summarized as follows. 
Including LL terms we get a 0.5% correction; the correction increases to 4% when higher-order 
terms proportional to 7^ are included; when all NLL contributions are included the correction 
decreases to 2.6%; when the Coulomb term is added the correction decreases to 2.1%; when the 
first NNLL terms (the Cjj constants) are included the correction grows to 4% for the choice A = 2, 
to 6% for the choice A = 0. As we have seen, the A = choice is somewhat disfavoured, since it 
overestimates the O(cv^) cross-section above threshold. We emphasize that the overall compression 
of the uncertainty band toward the high cross-section values is significant and always present once 
we include the NLL terms, regardless of the presence of the CV, and Coulomb terms, and regardless 
of our choice of the value of A. 

In the case of top productio at the Tevatron, it was shown in ref. J7[ that inclusion of the 
sole NLL terms proportional to 7^ provided an improvement in the agreement with the O(a^) 
partonic cross-section for the choice of scale \x = m t . This improvement is however accidental, 
and indeed, as shown in this paper, it is not sufficient to improve the accuracy of the resummed 



25 



calculation when the renormalization scale is changed to different values. Only at the full NLL 
level a significant reduction in the scale dependence is found. Similar results are obtained at higher 
energies, including the energy of the future LHC collider. 

The improvement in the predictive power allowed by the calculation we presented in this 
work is even more impressive when one considers the case of the bottom-quark production near 
threshold, at the future HERAB experiment. In this case the size of the strong coupling constant 
gives rise to very large perturbative corrections, which result in a factor of 4 uncertainty in the 
NLO result at a fixed value of the bottom mass and parton densities. Inclusion of resummation 
effects stabilizes the predicted cross-section at the level of ±50%. As in the case of the top 
cross-sections, the evaluation of the fully resummed cross-sections shows that the estimates of the 
theoretical uncertainty at the NLO obtained by varying the renormalization scales in the range 
mg/2 < fi < 2rriQ were correct. While the resummed result has a much smaller uncertainty than 
the fixed-order one, the resulting range of predictions is always included in the uncertainty band 
estimated at the NLO. 

Acknowledgement: We thank Matteo Cacciari for pointing out an error in Fig. 14 of the first 
version of this paper. 
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